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ABSTRACT 



Aims. We test predictions of evolution models on mixing of CNO-cycled products in massive stars from a fundamental perspective. 
Relative changes within the theoretical C:N:0 abundance ratios and the buildup of helium are compared with observational results. 
Methods. A sample of well-studied Galactic massive stars is presented. High-quality optical spectra are carefully analysed using 
improved NLTE line-formation and comprehensive analysis strategies. The results are put in the context of the existing literature data. 
Results. A tight trend in the observed N/C vs. N/0 ratios and the buildup of helium is found from the self-consistent analysis of 
main-sequence to supergiant stars for the first time. The catalytic nature of the CNO-cycles is confirmed quantitatively, though further 
investigations are required to derive a fully consistent picture. Our observational results support the case of strong mixing, as predicted 
e.g. by evolution models that consider magnetic fields or by models that have gone through the first dredge-up in the case of many 
supergiants. 
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1. Introduction 

Energy production in massive stars is governed by the CNO- 
cycles throughout most of their lifetimes. The general correct- 
ness of our understan ding of the CN O cycles and of the rele- 
vant nuclear data (e.g. iMaedeJflQSSl) is confirmed impressively 
by observation: when massive stars enter the Wolf-Rayet phase 
as WN subtypes, equilibrium CN O-processed ma terial becomes 
exposed on their surface (see e.g. ICrowtheill2007h . 

However, traces of mixing of CNO-cycled products from the 
stellar core to the stellar surface can already be found much ear- 
lier in the lives of massive stars. Observational indications of 
superficial abundance anomalies for carbon, nitrogen, and oxy- 
gen (and the burning product helium) in OB-type stars even on 
the main sequence (MS) and, more prominently, in the blue su- 
pergiants was found early in classification spectrograms (e.g. 
IWalbornlll976 ). Subsequent analyses provided evidence for a 
characteristic enrichment of nitrogen - which is the easiest to 
be detected - and helium in many e arly-type stars, both near 
the MS and in blue supergiants (e.g. 'Schonb erner et al Il988t 
iGies & Lambert 1992; Herrero et al. 1992; Ki lianlll99a IVennI 
ri995LILvubimkovll996HMcErlean et al.l ll999). 

A theoretical understanding of early mixing of CNO-pro- 
cessed material to the stellar surface could not be achieved 
within the framework of evolution calculations for non-rotating 
stars with mass-loss, which were state-of-the-art at that time (e.g. 
[Chio si & Maeder 1986). The pollution of the surface layers with 
CNO-cycled material only occurs when the star reaches the red 
supergiant phase in such models, via convective dredge-up. 
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* Based on observations obtained at the European Southern Obser- 
vatory proposals 62.H-0176 & 074.B-0455(A). 

** Based on observations collected at the Centro Astronomico Hispa- 
no Aleman at Calar Alto, proposals H2001-2.2-01 1 & H2005-2.2-016. 



Considered a secondary effect for a long time, stellar rotation 
has come lately into focus. It turned out that rotationally-induced 
mixing through meridional circulation and turbulent diffusion in 
rotating stars provides the means to change all model outputs 
substantially and to bring theory and observation into much bet- 
ter agreement (Maeder & Mevnet 2000; Heger & Langer 200(|). 
The latest step taken in the modelling was to consider the effects 
from an interplay of rotation and magnetic fields, which - de- 
pending on the detailed input physics and approximations made 
- only pr ovide minor modifications to the surface abundances 
(iHeger et a l. 2005) or substantial changes (Maeder & Meynej 



I2005L henceforth abbreviated as MM05). These differences re- 
sult mainly from the two groups using different sets of equations 
for computing the effects of magnetic fields; see in part icular, 
the changes brought to the system of equations given by I Spruit! 
(2002) by MM05 in their Sect. 2. 

The only means to verify the models is via systematic 
comparison with observations covering the relevant parts of 
the Hertzsprung-Russell diagram. Homogeneous analyses of a 
larger star sample from the main sequence to the supergiant stage 
have only recently become available. The results in particular 
for N abundances appar ently challenge the concept of rotational 
mixing in massive stars (Hunter et al.'2009) and thus the present- 
day evolution models (see, however, Maeder et al. 2009) . 

Here, we address the topic of early mixing of CNO-cycled 
products in massive stars from a fundamental perspective. We 
compare the model predictions for the relative changes within 
the C:N:0 abundance ratios and the buildup of helium (Sect. ^ 
with observations. For this we briefly review the status of the 
literature on CNO abundances and introduce a well-studied 
Galactic star sample for which high-accuracy analyses of high- 
quality spectra were performed using improved NLTE modelling 
and comprehensive analysis techniques (Sect. [3]l. Conclusions 
from this comparison are drawn in Sect. |4] 
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Fig. 1. Status of observational constraints on the (magneto-)hydrodynamic mixing of CNO-burning produ cts in massive stars from 
previous NLTE analyse s. Mass ratio s N/C over N/0 are displayed. Left pane l: main -sequen ce stars. Circles: Kilianl(ll992h: triangles: 
Gies& Lambert' ('1992'); diamonds: 'Cunha & Lambert' ('1994'), 'Daflon et al.' ('1999', '2001a"b^; squares: 'Morel et al.^ ('2008'); crosses: 
Hunter et al. (2 09). R ight panel: B A supergiant s. Triangles: Venn ( 1995), Venn & Przybilla (2003); circles: Takeda (2000); squares: 
Crowther et alj (l2006l) ; diamonds: ISearle etld] (1200 8). Error bars are omitted for clarity: uncertainties in the abundances of the 
individual elements are typically about a factor 2, such that the error bars can be lar ger than the plotting rang e. The lines represent 
predictions from evolution calculations, for a rotating 15 M© star (v™ - 300 km s~' . iMevnet & Maederll2003l MM03; until the end 
of the MS: solid red line, until the end of He burning: dashed blue line) and f or a star of the same mass and v™ that in addition takes 
the interaction of rotation and a magnetic dynamo into account (Maeder & Mevnet 2005. MM05; until the end of the MS: dotted 
line), respectively. The predicted trends are similar for the entire mass range under investigation, see Fig|3] 



2. Theoretical considerations 

Diagnostic plots of an abundance ratio vs. another one, such as 
N/C vs. N/0 or helium vs. N/0 discussed in this paper, depend 
on both the changes produced by nuclear reactions and the dilu- 
tion effects produced by mixing. Let us estimate the slope gk^jT 
produced by the nuclear effects at the beginning of CNO burn- 
ing, with everything expressed in mass fractions. Apart from the 
very massive stars (M > 40 Mq), one may assume that, at the be- 
ginning of the burning, the '^N enhancement comes from the '^C 
destruction via the CN cycle, and the oxygen content O remains 
about constant. Thus, one has dC - -^dN since ''*N globally 
results from the addition of two protons to '^C, 



d(N/0) = dN/0 



„ ,_ dN N dC , dN 
d(A^/C)=--^-dA.= - 



, 6 A^ 
'^7C 



(1) 



(2) 



This gives the slope 

d(N/C) ^ (N/C) 
d(N/0) {N/0) 

always in mass fractions. With initiaQ ratios N/C and N /O of 
0.31 and 0.11, respectively, we get 



, 6 A^ 
'^IC 



d(JV/C) 
d{N/0) 



3.11 



(3) 



This ratio is evidently greater than 1, since as A^ starts growing, 
C decreases, while O does not vary much. The relation turns 
slightly upward as N/C is increasing owing to the term in brack- 
ets in Eq. (|2]|. However, at some advanced stage in evolution. 



' The chemical composition of the evolution models is identi- 
cal to those used in the OPAL opacity tables: the solar mixture of 



iGrevesse & Noelsl ( Il993h . 



corresponding to WN stars not shown here, the curve will sat- 
urate and turn down slightly (Maeder 2009, p. 699), since the 
CN cycle is then at equilibrium, while '^O is still turned to '"*N. 
Dilution mixes a fraction f of N + AN enriched and C depleted 
materials with a fraction (1 - /) of the original A^ and C. Under 
the same assumptions as above, it is easy to show that, to the first 
order, the slope for the relative enrichments in the N/C vs. N/0 
plot behaves the same way as in Eq. (|2]l independently of /. 

The value of /, however, determines the amplitudes of the 
departures from the cosmic ratios. Our models with rotational 
mixing (Meynet & Maeder 2003, MM03; Ekstrom et al. 2008, 
E08) or with rotation and magnetic fields (MM05), as illus- 
trated e.g. in Fig. [1] have an initial slope J,J,2 ~4, which is 
in excellent agreement with Eq. ^. The amplitude / of the 
mixing depends on the various model assumptions, in particu- 
lar on the treatment of the shear mixing with or without hor- 
izontal tu rbulence. The m odels without horizontal turbulence 
(iMevnet & Maeder| l2000l MMOO) predict more mixing than 
models that account for it (MM03). Models that include both 
rotation and magnetic field predict a still larger mixing (MM05). 

Let us now consider the behaviour of the helium surface con- 
tent Ys vs. N/0 (as illustrated later). Strictly and only at the very 
beginning of the CN burning, and under the assumption of an 
initially constant oxygen, we get dFs = f dA^, since when 4 units 
of mass of helium are made, 14 units of mass of nitrogen are 
produced. The slope is 

-^II- ^-o^ 0.286 X 0.009 = 0.0026 , (4) 

d{N/0) 1 

i.e., it is essentially flat initially. Later in the evolution, both N 
and O change simultaneously, and one has to rely on numeri- 
cal models. The resulting slope in the models can vary; e.g., a 
steeper slope is inferred for the 25 Mq model of both MM03 
and models by MMOO than for the 15 Mq models of MM03 and 
MM05. This depends on whether the matter that arrives at the 



N. Przybilla et al.: Mixing of CNO-cycled matter in massive stars 



surface comes from inner regions that are at both CN and ON 
equihbria, or only at CN equilibrium. There are, of course, a 
range of intermediate cases. 

3. Observational constraints 

Numerous studies of CNO abundances in massive stars of the 
Milky Way are available from the literature, mostly for early 
B-type stars close to the MS and for B and A-type super- 
giants. We illustrate the results of several key publications in 
the N/O-N/C diagrams of Fig. [T] Sorne of the more recent stud - 
ies (ICrowther et alj|2006l; ISearle et al.ll2008l; iHunter et al.ll2009l) 
are based on NLTE model atmospheres, while the bulk of the 
data were obtained from NLTE line-formation computations on 
LTE model atmospheres - which is equivalent to the full NLTE 
appr oach in the cases under consideration (Nieva & Przvbilli 
I2OO7 . NP07). 

The MS stars show a concentration around the solar ratios 
of NIC ^ 0.3 and NjO a; 0. 1 . However, overall a wide range of 
N/O-N/C combinations has been realised, with the deviations 
from the predictions increasing in the supergiants. On the one 
hand, this leaves room for broad interpretation - even more so if 
only one of the elements (like N) is considered - including the 
statement that some observational data points pose a challenge 
for the evolution models. On the other hand, most of the data are 
nevertheless consistent with the predictions, as the abundance 
uncertainties are very large. Typically, the statistical Icr-error 
in abundance per element is about a factor ~2, and systematic 
uncertainties are often largely underestimated (for a discussion 
of this see Nieva & Przvbillaj 120101) or even unaccounted for. 
The error bars in Fig. [1] are larger than the entire plotting range 
in many cases. In consequence, no definite conclusions can be 
drawn on the quality of the stellar evolution models from these 
data. 

In past years a number of high-accuracy studies of mas- 
sive stars in the solar neighbourhood have been published 
by us. Carefully analysed data are availabl e on 6 slowly- 
rotating early B-type stars near the MS ( Nieva & Przvbillal 
[2OO6, 2007, 2008; Pr zybilla et al. 2008, PNB08) and on 14 
BA-type supergiants dPrzvbilla etaP l2006t iFirnsteinI l2006t 
ISchiUer & PrzvbiUall20(m 

In brief, high-resolution (R - 40-48 000) and high-S/N spec- 
tra (S /N>300) with wide wavelength coverage and thorough 
continuum normalisation (obtained with FOCES@Calar Alto 
Observatory and FEROS@ESO/La Silla) were analysed using 
a hybrid NLTE approach (Przybilla et al. 2006; NP07). State-of- 
the-art atomic input data were used in the modelling. In contrast 
to all previous work, multiple hydrogen lines, the helium lines, 
multiple metal ionisation equilibria and the stellar energy distri- 
butions were reproduced simultaneously in an iterative approach 
to determine the stellar atmospheric parameters. Chemical abun- 
dances were derived from analysis of practically the entire ob- 
servable spectrum per element. The rewards of such a com- 
prehensive, but time-consuming procedure are unprecedentedly 
small statistical error margins and largely reduced systematics. 
The relevant results - effective temperature Tgff, surface gravity 
log^, projected rotational velocity vsin/, surface helium abun- 
dance Y5, CNO abundances and their total, and the zero-age MS 
mass (derived under the assumption that the objects have evolved 
directly from the MS) - are summarised in Table [1] and visu- 
alised in Figs. I2H6] We re-ran the analysis of the supergiant tar- 
gets for the present work, taking advantage of improved model 
grids. Our new results for these supergiants agree with the earlier 
ones within the uncertainties but are more accurate. 
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Fig. 2. Sample stars in the log Teff-log g-plane. Diamonds denote 
MS stars and circles supergiants. The symbol size encodes the 
N/C ratio (by mass), with typical values exemplified in the leg- 
end. Evolutionary tracks for rotating stars of metallicity Z = 0.02 
(MM03) are shown. Predicted N/C ratios are indicated along the 
tracks. 



The sample is displayed in the log rgft-log ?-plane and com- 
pared to evolutionary tracks for rotating star^ in Fig. |2] The 
stars have initial masses between about 9 to 25 Mq. Stars close 
to the MS are apparently slow rotators, while the atmospheric 
expansion has erased any indication of the initial rotational ve- 
locity in the supergiants. There is good qualitative agreement be- 
tween the observations and predictions for the N/C ratios, find- 
ing low values close to the initial (solar) ratio near the MS and 
enhanced values in the supergiants. However, the observed N/C 
ratios reach much higher values than predicted by the models 
discussed in Fig. |2] provided these stars have recently left the 
MS and evolve now into red supergiants. 

The behaviour of light element abundances of the MS sam- 
ple stars is shown in Fig. |3] Theoretical evolutionary tracks are 
overplotted. It is striking to see that, regardless of the initial mass 
and of the initial rotation, the slope of the theoretical tracks is al- 
ways the same and equal to the nuclear effect of the CN cycle 
as derived analytically in Sect. |^ Thus the observed slope is a 
confirmation of the activity of the CN cycle. It does not depend 
on the stellar model but constitutes observational evidence of the 
nuclear path. As already indicated before, what is model depen- 
dent are the amplitudes of the departure from the cosmic ratios. 
Figure|3]shows that the observed points may be reproduced from 
models having differing velocities and masses, with and without 
magnetic field; however, by itself this diagram cannot help in 
disentangling the various possibilities. 



- The metal content (by mass) of the models is Z = 0.020, while a 
value of 0.014 seems more appropriate for the massive star population 
in the solar neighbourhood (PNB08). No fundamental changes to our 
conclusions are expected to arise from this difference. The tracks will 
shift a bit, in particular the zero-age MS towards higher gravities. 

^ The theoretical slope varies between ~3.0 and 4 .2 for typical ini- 
tial CNO abundances (e.g. solar values according to Ander s & Grevessg 
1989; Grevesse & Sauval 1998; Asplund et al. 2009, or the cosmic 
abundance standard of PNB08), indicated by the grey area in Fig. [3] 
Improved agreement between observation and theory may be achieved 
when tailored abundances are used for the stellar evolution models, see 
also the previous footnote. 
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Table 1. Stellar parameters and elemental abundances of the sample stars 



Star 


7-eff 

(K) 


(cgs) 


V sin ;■ 
(kms-') 


Ks 
(by mass) 


E{cy 


£(N) 


£(0) 


£(i;CNO) 


M/Me 


B-type main 


-sequence stars 


















HD36591 


27000+300 


4.12+0.05 


12+1 


0.28+0.03 


8.33+0.08 


7.75+0.09 


8.75+0.11 


8.92+0.07 


13 


HD61068 


26300+300 


4.15+0.05 


14+2 


0.28+0.03 


8.27+0.07 


8.00+0.12 


8.75+0.09 


8.93+0.06 


12 


HD63922 


31200+300 


3.95+0.05 


29+4 


0.25+0.03 


8.34+0.08 


7.77+0.08 


8.79+0.10 


8.95+0.07 


19 


HD74575 


22900+300 


3.60+0.05 


11+2 


0.28+0.03 


8.37+0.10 


7.92+0.10 


8.80+0.08 


8.98+0.06 


12 


HD 122980 


20800+300 


4.22+0.05 


18+1 


0.28+0.03 


8.32+0.09 


7.76+0.08 


8.72+0.05 


8.90+0.04 


8 


HD 14943 8 


32000+300 


4.30+0.05 


4+1 


0.28+0.03 


8.30+0.12 


8.16+0.12 


8.77+0.08 


8.97+0.06 


17 


B A- type supergiants: 


















HD 13476 


8500+150 


1.40+0.10 


12+2 


0.36+0.08 


8.18+0.11 


8.58+0.04 


8.63+0.06 


8.98+0.03 


16 


HD14433 


9150+150 


1.40+0.10 


17+3 


0.32+0.07 


8.23+0.04 


8.23+0.03 


8.67+0.05 


8.91+0.03 


19 


HD34085 


12100+150 


1.75+0.10 


25+3 


0.32+0.04 


8.23+0.09 


8.47+0.06 


8.75+0.05 


9.01+0.03 


23 


HD87737 


9600+150 


2.00+0.10 


2+2 


0.37+0.03 


8.25+0.06 


8.52+0.08 


8.73+0.06 


9.02+0.04 


10 


HD92207 


9500+200 


1.20+0.10 


30+5 


0.35+0.06 


8.33+0.08 


8.25+0.04 


8.79+0.07 


9.00+0.04 


27 


HD111613 


9150+150 


1.45+0.10 


17+2 


0.35+0.05 


8.29+0.10 


8.46+0.04 


8.73+0.04 


9.01+0.03 


17 


HD 187983 


9300+250 


1.60+0.15 


15+6 


0.32+0.05 


8.29+0.09 


8.19+0.04 


8.78+0.02 


8.98+0.02 


15 


HD 195324 


9200+150 


1.85+0.10 


3+3 


0.38+0.06 


8.10+0.11 


8.70+0.08 


8.73+0.04 


9.07+0.04 


11 


HD 197345 


8700+150 


1.20+0.10 


10+9 


0.37+0.04 


8.09+0.07 


8.56+0.07 


8.69+0.04 


8.99+0.03 


22 


HD202850 


10800+200 


1.85+0.10 


14+5 


0.38+0.06 


8.16+0.04 


8.70+0.06 


8.76+0.05 


9.09+0.03 


15 


HD207673 


9250+100 


1.80+0.10 


1+2 


0.33+0.08 


8.17+0.09 


8.48+0.03 


8.72+0.04 


8.99+0.03 


12 


HD2 10221 


8400+150 


1.40+0.10 


0+0 


0.34+0.03 


8.22+0.06 


8.52+0.06 


8.70+0.05 


9.00+0.03 


15 


HD212593 


11200+200 


2.10+0.10 


6+2 


0.35+0.05 


8.30+0.08 


8.44+0.06 


8.74+0.04 


9.01+0.03 


12 


HD2 13470 


8400+150 


1.30+0.10 


13+1 


0.32+0.07 


8.17+0.08 


8.54+0.04 


8.65+0.03 


8.97+0.02 


18 
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Fig. 3. Comparison of N/C vs. N/O abundance ratios (by mass) 
for our six MS sample stars (diamonds) with model predictions 
for 9 Mq (left) and 20 Mq stars (right panel) at different v™ (see 
legend, equivalent to 30, 50 and 70% of the star's breakup ve- 
locity, EOS). The dotted line (right panel only) describes the 
magnetic ISMq model of MM05 with v™ =300kms-'. The 
symbol size encodes the stellar mass and error bars give Icr- 
uncertainties. The long-dashed lines correspond to the theo- 
retical slope 3.77 as deduced in Sect. |2l The grey area spans 
the full range of theoretical slopes using different references 
fo r solar abundances (s ee text, lower envelope defined by data 
of lAsplund et al J 120091) and the cosmic abundance standard of 
PNB08 (upper envelope). 



The same evolutionary tracks are plotted in the N/C\s. log g 
diagram in Fig.|4] We see that, except for the stars HD 61068 and 
HD 149438 (t Sco, which is a genuine slow rotator with a rather 




Fig. 4. N/C abundance ratios on the MS as function of log^. 
The same tracks and observations as in Fig. |3] (right panel) are 
displayed. 



strong magnetic fielqj of probably fossil origin, iDonati et al.l 
2006), the other points may be accounted for by models having 
low initial rotation. HD 61068 is marginally compatible with the 
15 Mq model at v™ - 300 km s"' computed with magnetic field, 
while T Sco challenges present stellar models. It has a behaviour 
that may be explained by a homogeneous evolution, but this still 
has to be confirmed by further computations. This is a very inter- 
esting star that certainly deserves further inspections both from 
the observational and theoretical points of view. 

The behaviour of light element abundances in the whole star 
sample is shown in Fig. |5] In contrast to the literature values 
(Fig.[D> a clear and tight trend is found, confirming the predicted 
locus of N/O-N/C abundance ratios. However, as already indi- 



A magnetic field is present also in HD 74575 dHubrig et alj200jt) . 
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Fig. 5. N/C vs. N/0 abundance ratios (by mass) for our sam- 
ple stars. B-type MS stars are displayed as diamonds, BA-type 
supergiants as circles. The symbol size encodes the stellar mass 
and error bars give Icr-uncertainties. The different lines describe 
model predictions for 15 M© stars identical to those in Fig.[T] 



cated above, the MM03 models for rotating stars with mass loss 
evolving towards the red supergiant stage (solid line in Fig. |5]l 
predict mixing that is too low, i.e. too low / (Sect. |2]i, in partic- 
ular for most of the supergiants. A combination of five reasons 
may provide an explanation. 

i) Higher than average rotation velocities in the progenitor stars 
of these supergiants on the MS may reconcile the situation for 
some objects. 

ii) Evolution models for rotating stars that also account for the 
interaction of rotation and a magnetic dynamo (MM05) predict 
enhanced mixing signatures of the amount required (dotted line). 
Ill) Some stars may have evolved in a close binary, which can 
also lead to enhanced mixing associated with mass transfer, 
iv) Some objects may have been siblings to rSco on the MS, 
climbing up the N/O-N/C relation even further in their further 
evolution. 

v) Supergiants may already have evolved through the red super- 
giant phase (e.g., on a blue loop) to expose first dredge-up abun- 
dance ratios, which could quantitatively also explain the obser- 
vations (dashed line). 

More information may be derived from the helium content, 
which in the case of BA-type supergiants is determined here for 
a significant number of stars for the first time in a self-consistent 
analysis. Our results are displayed in Fig. |6] On the MS no he- 
lium surface enrichment is observed, as predicted in the models 
for stars with masses below about 20 M©. After the MS, the pic- 
ture is blurred by the possible occurrence of a blue loop. Actually 
the interpretation of the blue supergiant can become really con- 
straining only when we obtain additional hints to the previous 
evolution of the star. Has the blue supergiant evolved directly 
from the MS, or has it evolved in that stage after going through 
a red supergiant stage? At the moment, from the models the sit- 
uation would be the following: for models below 20 M©, He- 
enrichments at the level observed in the present supergiants are 
only compatible with models having undergone a dredge-up in 
the red supergiant phase. This is true whether rotation is con- 
sidered or not, or a magnetic field is accounted for or not. The 
present track for the magnetic 15 Mq model was computed only 
up to the end of the MS phase and thus did not yet go through the 
dredge-up phase. Depending on the rotation velocity, the pres- 
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Fig. 6. Surface helium abundance Fs as a function of N/0 ratio 
(both by mass) for our sample stars. Symbol and line encoding 
as in Fig.[T] 



ence of a magnetic field or its absence, models, after the red 
supergiant phase, will populate diverse parts of the region in the 
plane Y^ versus N/0, as illustrated e.g. by the dashed and dotted 
lines in Fig.|6] 

We note, however, that we cannot exclude at present the pos- 
sibility that the observed heUum abundances in the supergiants 
may be overestimated. A systematic downward shift by a mere 
10-15% (which is within the typical systematic uncertainties in 
our abundance determinations) would be sufficient, e.g., to bring 
the observations and the magnetic model in Fig. |6]into agree- 
ment. All observed lines arise from two energetically close He i 
levels only (2p'P°,^P°). Model atom shortcomings (such as in- 
sufficient ab-initio collisional data), which may become impor- 
tant only in the weak-line limit at the cool temperature border, 
could therefore remain unnoticed and could give rise to system- 
atics. It would not be the first time that uncertainties in atomic 
data complicate the statistic equilibrium and radiative transfer 
calculations neede d to interpret th e observed Hei lines (e.g. 
Przybilla 2005; Naiarro et al.l 120061). Further investigations are 
required before firm conclusions on the evolutionary state of the 
supergiants are drawn from the helium abundances, and a fully 
coherent picture can be deduced. 

Finally, an inspection of the total of the CNO abundances 
(see Table [U is instructive. Because the CNO cycles are cataly- 
tic, the sum of particles is conserved. Indeed, the sample 
shows only a small scatter around the mean combined abun- 
dance, log SCNO-H 12 = 8.98+0.05. This is further (indirect) 
evidence that chemical inhomogeneities among the massive 
star population in the solar neighbourhood are small (Nieva & 
Przybilla 2008; PNB08), even at extended distances out to 1- 
2 kpc from the Sun. Predictions of a high efficiency of hydro- 
dynamic mixing in the interstellar medium and consequently of 
a loca lly well-mixed Milky Way (iEdmundslll975HRov & Kunflil 
Il995h are thus supported as well. 

4. Conclusions 

Evolution models for massive stars predict very tight relations 
for the change of surface N/0 and N/C abundance ratios and for 
the buildup of helium, as a consequence of mixing with CNO- 
cycled matter (for given initial chemical composition). Massive 
stars in the solar neighbourhood, which supposedly share a 
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nearly uniform initial composition (PNB08), are therefore ex- 
pected to follow the predicted relations, which are governed by 
nuclear reactions and the dilution effects produced by mixing. 

A comparison with NLTE spectral analyses from the litera- 
ture on CNO abundances in massive stars of ~8-25 M© leaves 
room for broad interpretation because of large uncertainties. 
Observations may even be interpreted as posing a challenge to 
theory, in particular when the full set of information is not ac- 
counted for and conclusions are drawn from one indicator alone, 
e.g. nitrogen. 

On the other hand, our high-precision analyses of a sample 
of Galactic massive stars from the main sequence to the super- 
giant stage find these tight correlations. Even though further in- 
vestigations are required to refine the observational constraints 
on helium abundances before a fully coherent picture can be ob- 
tained, the case for strong mixing is clearly supported. It is pre- 
dicted, e.g., by the models of MMOO, by models with magnetic 
field (MM05), or in cases where the models have gone through 
the first dredge-up. The separation of the different possibilities 
may come from studing the evolution of the rotation velocities 
with time. On the one hand, models with magnetic field predict 
higher rotational velocities at the end of the MS phase, because 
the strong internal coupling transmits some of the fast core rota- 
tion to the surface. On the other hand, stars on a blue loop should 
be slower rotators than objects on their first passage from blue to 
red on average because of the additional angular-momentum loss 
experienced through strong mass loss during the red supergiant 
phase. Indeed, some of the slowest rotators show very strong 
mixing signatures (Table [TJ. 

Finally, we have to point out that some of our results have the 
potential of challenging the currently available evolution mod- 
els. The star tSco (HD149438) stands out in the sample as it 
shows characteristics that may be explained by a homogeneous 
evolution, but it requires a highly-efficient spin-down mecha- 
nism. One may speculate on magnetic breaking due to angular- 
momentum losses by a magnetically confined line-driven stellar 
wind or magnetic coupling to the accretion disc during the star- 
forming process in the case of a fossil field. Even though the 
topic is not understood theoretically i n a comprehensive w ay, 
spin-dow n times of the order of 1 Myr dUd-Doula et al.ll2009h or 
even less dMikulasek et al.l2008l) are reported for some magnetic 
massive stars, possibly leading to the required slow rotation al- 
ready on the zero-age MS in this case as well. Then, if the present 
helium abundances are confirmed and the supergiants are shown 
to have evolved directly from the MS, a different kind of mixing 
may also be required. 

We conclude that the tight observational constraints that are 
required for a thorough testing of the stellar evolution models are 
within reach. Applications of the improved modelling and analy- 
sis techniques to high-quality observations of larger star samples 
in the Milky Way and other galaxies may finally provide the em- 
pirical basis to benchmark the models. It may thus become fea- 
sible to disentangle the effects of metallicity, rotation, magnetic 
fields, and binarity on massive star evolution. 
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